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Abstract 

A particle rotor model is developed which couples several valence protons 
and neutrons to a rigid triaxial rotor core. It is applied to investigating the 
chirality in odd-A nucleus ^^^Nd with T^h\^i2 ® configuration for the 

first time in a fully quantal approach. For the two chiral sister bands, the 
observed energies and the B{M1) and B{E2) values for the in-band as well 
as interband transitions are reproduced excellently. Root mean square values 
of the angular momentum components and their probability distributions are 
used for discussing in detail the chiral geometry of the aplanar rotation and 
its evolution with angular momentum. Chirality is found to change from a 
soft chiral vibration to nearly static chirality at spin / = 39/2 and back to 
another type of chiral vibration at higher spin. 
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Spontaneous chiral symmetry breaking is a phenomenon of general inter- 
est in chemistry, biology and particle physics. Since the pioneering work of 
nuclear chirality in 1997 [l| , much effort has been devoted to further explore 
this interesting phenomenon. Following the observation of chiral doublet 
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bands in = 75 isotones [2j, more candidates have been reported over more 
than 20 nuclei experimentally in A~ 100, 130 and 190 mass regions including 
odd-odd, odd-A and even-even nuclei 0, i, i, S, 0, H i, 0, H Q 13 



Chiral symmetry breaking was initially suggested to occur in a stable tri- 
axial deformed nucleus, with a high-j particle-like valence proton (neutron) 
and a high-j hole-like valence neutron (proton) [1] . In this case, the angular 
momenta of the core, the valence proton and neutron are mutually perpen- 
dicular, arranging in the body-fixed frame into two systems with opposite 
chirality (left- and right-handed). Due to the restoration of chiral symmetry 
by quantum tunneling, a pair of near degenerate AI = 1 bands -the chiral 
sister or doublet bands- are observed [H, IHj . 

Theory wise, chiral doublet bands were first investigated in the one- 
particle-one-hole-rotor model (PRM) and the corresponding tilted axis crank- 
ing (TAG) approximation [ij. Later on realistic TAG approaches, as the 
Strutinsky shell correction method with a hybrid Woods-Saxon and Nilsson 
potential [l6|, the Skyrme Hartree-Fock model 17|, as well as the relativistic 



mean field model |18l . |19| has been developed to investigate this new phe- 
nomena. Within the TAG mean field approximation, the left-handed and 
right-handed solutions are exactly degenerate. It is not possible to calcu- 
late the energy difference between the bands, which is the consequence of 
quantum tunneling between the two solutions. Before the onset of chirality 
of the mean field, the precursor of the symmetry breaking occurs as a soft 
vibration between the right- and left-handed configurations. These chiral 
vibrations have been studied in the framework of the random phase approxi- 
mation (RPA) based on the TAG mean field 0,[2l|. However, TAG+RPA is 
not able to describe the smooth transition from a slow vibration to quantum 
tunneling between the left- and right-handed mean field solutions. Moreover, 
angular momentum is not a good quantum number and the electromagnetic 
transitions are calculated in semiclassical approximation. 

The fully quantal PRM does not suffer from these deficiencies. Total 
angular momentum is a good quantum number. The energies and transi- 
tion probabilities are treated fully quantal. In particular, it provides the 
energy splitting between doublet bands covering the whole range from chiral 
vibrations to weak tunneling between the two chiral configurations. On the 
other hand, the PRM is based the assumption that the rotor has a fixed 
deformation. This may appear as a problem, because the nuclei in which 
chiral doublets have been found are considered to be soft with respect to the 
triaxiality parameter 7. The study in the framework of an IBA core cou- 
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pled to a particle and a hole found a substantial coupling between the shape 
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In contrast, the 
indicate that 



and angular momentum orientation degrees of freedom 
microscopic studies in the framework of TAC+RPA [2 
the angular momentum dynamics is almost decoupled from from the shape 
degrees of freedom, which justifies the application of the PRM. 

For odd-odd nuclei in A ~ 100 and 130 regions, chirality has been exten- 
sively studied with PRM with 1-particle-l-hole configuration 23,|23], or by 
introducing pairing to simulate the effect of many valence nucleons [gI, [2^, 
26l . |27| | . Good agreement with the experimental spectra and electromagnetic 



transitions has been obtained by PRM for the doublet bands in odd-odd 
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nuclei i26Cs ^^SQs gQ, and ^^^Rh 

Apart from odd-odd nuclei, chiral doublet bands have been observed in 
odd-A nuclei. The first example is the chiral spectra characteristic observed 
in ^^^Nd which is further confirmed by the lifetime measurements [2(1] and 
suggested to be built on the configuration vr/i^^^g ® ^^11/2 • More candidates 
have been reported in ^^^Rh [3, ^°^Rh 0,[l3 as well. For even-even nuclei, a 
candidate was suggested for ^^^Nd j^, however the interpretation is doubted 



on the basis of the recent lifetime measurement [29|. The existence of these 
data makes a PRM that treats more than one valence proton and one valence 
neutron highly desirable. A PRM with an axial or a triaxial core coupled to 
many-particle configurations has been developed and used for the description 
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of magnetic bands of i^S'i^^Pb 30] and the wobbling excitation in ^^^Lu 
However the study of chirality with such kind of model is still absent. There 
has been a dispute on the identification and interpretation of chiral doublet 



bands based on 1-particle-l-hole configuration in odd-odd nuclei j22|, |28|, l29|. 



32I . I33I ] . The study of the doublet bands in odd-A nuclei in the framework of 
a rotor coupled to two particles and one hole should shed new light on the 
question of nuclear chirality. 

In this Letter, a triaxial n-particle-n-hole PRM will be developed to treat 
more than one valence proton and one valence neutron and applied to the 
study of nuclear chirality. The energy spectra and electromagnetic tran- 
sitions probabihties of the doublet bands in ^^^Nd will be calculated and 
compared with the data available as well as the previous TAC-I-RPA results. 
The structure of doublet bands with 2-particle-l-hole configuration will be 
discussed. 

The total Hamiltonian is expressed as. 



H = ^coll + 



'intr' 



(1) 
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with the collective rotor Hamiltonian H^^^i^ 



^-11-1.^ = 1.^;^' (2) 

where the indices k = 1,2,3 refer to the three principal axes of the body-fixed 
frame, R^, Ik, Jk denote the angular momentum operators for the core, the 
total nucleus and the valence nucleons, respectively. The moments of inertia 
for irrotational flow are adopted, i.e., J'k = Jq sin^(7 — 2Tik/?i). The intrinsic 
Hamiltonian for valence nucleons is 

■^intr ~ ^ ^ ^p,u(ip^u^p,u + ^ ^^ ^nyo^ny^^y 1 (3) 

u v' 

where Sp^y and Sn^v' are the single particle energy of proton and neutron. 
The single particle states are expressed as 

«;r|0) = 5^cai«,fi), a+|0) = 5^(-ir^ci';i|a,-fi), (4) 

where f2 is the projection of the single-particle angular momentum j along 
the 3-axis and is restricted to half of the values, — j, ■ ■ ■ , — 1/2, 1/2, 3/2, ■ ■ ■ ,j, 
due to the time- reversal degeneracy [3J|, and a denotes the other quantum 
numbers. For a system with z valence protons and n valence neutrons, the 
intrinsic wave function is given as 

I*') = (ft "l) (ft <^ (ft -iy) [U 10) (5) 

with zi + Z2 = z,ni + n2 = n, < zi < z,0 < rii < n. 

The total wave function can be expanded into the strong coupling basis, 

\IM) = J2cKMMKip), (6) 

Kip 



with 



\IMK^) = {\IMKM + (-1)^-^|/M - K)m) , (7) 

^^2(1 + dicoV^) 
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where \IMK) denotes the Wigner functions y and is a shorthand 

notation for the configurations in Eq. ([5]). The basis states are symmetrized 
under the point group which leads to K — | (2:1 — Z2) — \ {ni — 722) being 
an even integer with Vt = ■ ■ ■ , —3/2, 1/2, 5/2, ■ ■ -. The reduced transition 
probabihties B{M1) and B{E2) can be obtained from the wave function of 
PRM with the Ml and E2 operators 

The single particle energy for proton and neutron Ep^i, and Eny in Eq. ([3]) 
are provided by the triaxial deformed single-j shell Hamiltonian [l[. 



ksp = ±\C [co.jGI - + ^(ii +m'>, (8) 

where the plus or minus sign refers to particle or hole, and the coefficient C 



is proportional to the quadrupole deformation (3 as in Ref. [271 . 

In the PRM calculations for the doublet bands in ^^^Nd, the configuration 
^^11/2 ® ^^11/2 l'^^] adopted. The deformation parameters P = 0.235 
and 7 = 22.4° for ^^^Nd are obtained from the microscopic self-consistent 



triaxial relativistic mean field calculation [18[. Accordingly, Cp and Cn in 
Eq. ([8]) take values of 0.323 and —0.323 MeV 27|. The moment of inertia 
Jo = 29.0 MeV/h'^ is adjusted to the experimental energy spectra. For 
the electromagnetic transition, the empirical intrinsic quadrupole moment 
Qo = (3/v^57r)i?gZ/5 = 4.0 eb, gyromagnetic ratios = Z/A = 0.44, and 
(7p=1.21, 5f„=-0.21 are adopted |25|. 

The calculated excitation energy spectra E{I) for the doublet bands A 
and B in ^^^Nd are presented in Fig. [H together with the corresponding 
data 0, [i^] ■ The experimental energy spectra are excellently reproduced by 
the PRM calculation. Except for the 27/2-state in band B, the calculated 
results agree with the data within 50 keV. In particular, the trend and am- 
plitude for the energy splitting between two partner bands are excellently 
reproduced. 

For comparison, the TAC+RPA results are also included. The energy 
splitting between two sister bands is well described for I < 39/2, which is the 
vibrational regime, where the RPA is stable. The energy splitting disappears 
at / = 41/2, where the TAG solution attains chirality. The description of 
the experiment in the region J > 39/2 is beyond the realm of RPA based on 
the planar TAG solution. Being a quantum theory that is not restricted to 
small amplitude vibrations, PRM is able to perfectly reproduce the energy 
splitting for the whole observed spin region, staying within the configuration 
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^^11/2 ® ^^11/2 triaxial rotor. The analysis below shows that the 

motion of the angular momentum vector in fact changes from chiral vibration 
about the plane spanned by the long and short axes of the triaxial shape to 
tunneling between the left- and right-handed configurations. 

The calculated in-band and interband transition probabilities B{M1) (up- 
per panel) and B{E2) (lower panel) for the doublet bands in ^^^Nd are pre- 
sented in Fig. [21 together with the available data and the TAC+RPA re- 
sults j2o|. For / < 39/2, the observed in-band B(M1) values in the two 
bands are almost the same and much larger than the interband ones. These 
features are perfectly reproduced by the PRM calculation. The TAC+RPA 
gives a similar result. After I = 39/2 the two models differ. In the case 
of PRM, the interband and in-band B{M1) transitions become comparable 
due to the transition to tunneling regime, which is absent in the TAC+RPA 
calculations. 

For / < 39/2, the observed in-band B{E2) values, in particular their sim- 
ilarity, are well reproduced by both the PRM and TAC+RPA calculations, 
while the observed interband B{E2) values are underestimated in both the 
PRM and TAC calculations. The kink at / = 39/2 of the in-band B{E2) 
values of the PRM reflects the transition to tunneling regime as mentioned 
in the discussion of B{M1). The ratio between the in-band and interband 
B{E2) and B{M1) values depends sensitively on the details of the transi- 
tion from the vibrational to the tunneling regime, which may account for the 
deviations of the PRM calculation from experiment. 

The success in reproducing the energy spectra and transition probabilities 
for the doublet bands A and B in ^^^Nd suggests that the PRM calculation 
must correctly account for the structure of the states. To exhibit their chiral 



geometry [26|, we calculated for the bands A and B in Nd the expectation 



values of the squared angular momentum components for the total nucleus 



Ik = \/(/|), the core Rk = \/ (Rl), the valence neutron Jnk = \ (ilk) 



and the valence protons Jpk = y {{j{pi)k + j(p2)fc)^)) {k = 1,2,3), which are 
presented in Fig. [3] and HI 

As shown in Fig. [3], for both bands A and B, the collective core angular 
momentum mainly aligns along the intermediate axis (i-axis), because it has 
the largest moment of inertia. For 7 = 22.4° the ratios between the moments 
of inertia are Ji : Js '■ Ji = 6.8 : 2.6 : 1.0. The angular momentum of the 
/iii/2 valence neutron hole mainly aligns along the long axis (/-axis) and that 
of the two /iii/2 valence protons mainly along the short axis (s-axis), which 
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correspond to the orientation preferred by their interaction with the triaxial 
core To be more precise, Jn ~ 5h along /-axis, and both Jp ~ lOh and 
R ~ 6 — 12h lie in the plane spanned by i- and s- axis, which together form 
the chiral geometry of aplanar rotation. 

In order to show the picture more clearly and examine the evolution of 
the chiral rotation, the root mean square values of the angular momentum 
components for the nucleus, the core, the valence neutron, and protons at spin 
I = 29/2, 39/2, and 45/2 are illustrated in Fig. HI The i- and s- components 
of Jn, and the /- components of Jp and R, which are negligible, have been 
ignored for clarity. As the total angular momentum increases, R increases 
gradually, Jn remains almost unchanged, while Jp moves gradually toward 
the z-axis. The difference between proton and neutron is due to alignment 
effect by the Coriolis force, which is much weaker for the single neutron hole. 
At spin / = 39/2, where the doublet bands have smallest energy difference, 
the orientations of R, Jn and Jp for band A and band B are nearly identical. 
Hence around this spin the structure comes closest to the ideal chiral picture 
of a left- and a right-handed configurations with equal components R, Jn, 
and Jp along the respective i-, /-, and s-axes, which weakly communicate. 

Further insight into the development of chirality with increasing angular 
momentum provides Fig. [5l which shows the probability distributions for the 
projection of the total angular momentum along the i- and s-axes. For 
the PRM state Eq. the probability for the projection K of total angular 
momentum on quantization axis is 



For triaxiality parameter 7 = 22.4°, the /-axis is used for quantization. The 
distributions with respect to the other axes are obtained by placing 7 into 
another sector, such that the shape is the same but the principal axes are 
exchanged. The probability distribution of the angular momentum projection 
on the 2-axis is given by Eq. (Q and 7 = 97.6°; and the distribution of the 
projection on the s-axis is obtained with 142.4°. 

For spin I = 29/2, near the band head, the probability distribution of 
two bands differ as expected for a chiral vibration. For the lower band A the 
maximum probability for the i-axis appears at Ki = 0, whereas the proba- 
bility for the higher band B is zero aX Ki = 0, having its peak at Ki = 25/2. 
These are the typical probability distributions expected for the zero phonon 
state (A), which has a wave function that is symmetric in the projection on 




(9) 
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the i-axis, and the one phonon state (B), which is antisymmetric in the same 
degree of freedom. The probabihty distributions with respect to the /-axis 
have a peak near Ki = 11/2, which is generated by the /111/2 neutron hole. 
The probabihty distributions with respect to the s-axis have a peak near 
Kg = 8, which is generated by the pair of /iii/2 protons. Hence the chiral 
vibration consists in an oscillation of the collective angular momentum vector 
R through the s-/-plane. This reveals the structure of the chiral vibration 



calculated by means of TAC+RPA in Ref. [20 . 

At spin I = 39/2, the probability distributions for band A and B are very 
similar. The distributions are peaked at Ki = 7, Ki = 17, Ks = 17, which 
corresponds to a vector with an orientation close to the vector I in Fig.m The 
distributions show the characteristics of static chirality. The finite values of 
P{Ki = 0) and P{Ks = 0) reflect the tunneling between the left- and right- 
handed configuration, which is responsible for the remaining 76 keV energy 
difference between the bands. The tunneling is two-dimensional: along the i- 
and s-axes. If it was only along the i-axis, band B should have P{Ki = 0) = 0, 
as for the chiral vibration. The well developed tunneling regime is restricted 
to / = 39/2. For higher spin, where the energy difference between the chiral 
partners increases, they attain vibration character again. This is reflected 
by the increasing differences between the probability distributions of bands 
A and B. In particular the fact that P{Ks = 0) is finite for band A (zero 
phonon) and zero for band B (one phonon) shows that the motion contains a 
vibration of the vector I through the /-i-plane. The two vibrational regimes of 
our quantal calculation agree with the ones found by means of the TAC+RPA 
calculations in Ref. (2l| . 

In summary, we have developed a particle-rotor model, which couples a 
triaxial rotor with many valence protons and many valence neutrons. We 
used it to investigate the chirality of the configuration T^hf^^^ ® ^^11/2 
odd-A nucleus ^^^Nd. The energy spectra of the doublet bands, the reduced 
probabilities B{M1) and B{E2) for in-band as well as for the interband tran- 
sitions are reproduced excellently. Remarkably, the agreement with the data, 
in particular the energy splitting, is considerably better than in the case of 



odd-odd nuclei [23[, indicating that the longer angular momentum vector of 
the two protons leads to more pronounced chirality. The chirality is shown 
to be a transient phenomenon. The chiral partner bands start as a soft vi- 
bration of the angular momentum perpendicular to the plane spanned by the 
short and long axes, where band A is the zero- and band B the one-phonon 
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state. With increasing angular momentum the vibration becomes strongly 
anharmonic, progressively localizing in left- and right-handed configurations. 
Maximal chirality is reached at / = 39/2, where the two bands approach each 
other closest. At this spin they have very similar distributions of the angu- 
lar momenta of the core, the valence neutron, and the valence protons in 
the left- and right-handed sectors, which reflect reduced left-right tunneling. 
With further increasing spin the two bands again develop into the zero and 
one-phonon states of a chiral vibration of the angular momentum about the 
intermediate axis. The success of the present model encourages further ap- 
plication for chiral rotation in other nuclei including the even-even nuclei and 
settling the dispute on the identification and interpretation of chiral doublet 
bands based on PRM with 1-particle-l-hole configuration. 
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Figure 1: (color online) The excitation energies E{I) for the chiral sister bands in ^^^Nd 
calculated by means of the triaxial PRM with configuration nh'^^^^'^i/h^^^^ (open symbols) 
in comparison with the data (filled symbols) 0, HO] and the corresponding TAC+RPA 
results (dotted lines) [2^. The energies are relative to the band head Eq of the chiral 
bands, with a rotor reference subtracted. 
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Figure 2: (color online) The B{M1) and B{E2) values calculated by means of the PRM 
for the chiral sister bands in ^^^Nd (open symbols) in comparison with the data (filled 
symbols) [20] and the corresponding TAC+RPA results (dotted lines) (20| . 
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Figure 3: (color online) The root mean square components along tlie intermediate (i-, 
circles), short (s-, squares) and long {1-, triangles) axis of the core = \J {R\) , valence 

neutron J^k = \l (j^fc), and valence protons angular momenta Jpk = \j {{j{pi)k + i{p2)kY) 
calculated as functions of spin / by means of the PRM for the doublet bands in ^^^Nd. 
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Figure 4: (color online) The angular momentum geometry in PRM for the doublet bands 
in ^35Nd at spins / =29/2, 39/2, and 45/2. 
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Figure 5: (color online) The probability distributions for projection of total angular 
momentum on the long (/-), intermediate (i-) and short (s-) axis in PRM for the doublet 
bands in ^^^Nd. 



16 



